Localization of electronic states by fullerene charges in carbon nanotubes 
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We study the effects of the electrostatic interaction produced by charged fullerenes encapsulated 
in carbon nanotubes, showing that they are able to modify locally the electronic density of states 
in the hybrid system. In the cases where the interaction is felt as an attractive potential by the 
electrons in the nanotube, localized electronic states are formed in the nanotubes around the position 
of the fullerenes. This produces an effective narrowing of the gap in semiconducting nanotubes over 
a distance of a few nanometers, in agreement with the spatial modulation of the gap observed in 
the experiments. 
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Carbon nanotubes have nowadays a great potential for 
technological applications in devices at the nanometer 
scale. This makes very important the precise knowl- 
edge of their electronic properties, once the first step 
of discerning their metallic or insulating behavior has 
been completed after confirmation of the theoretical 
analyses pj by experimental results 0- We have learned 
from these studies that the low-energy spectrum of a car- 
bon nanotube may be gapless or not, depending on the 
helicity of the tubular structure. 

Hybrid structures have been recently synthesized by 
encapsulation of fullerenes in the hollow space of the 
nanotubes. It has been shown that the insertion of the 
fullerene molecules in these so-called peapods leads to a 
significant modification of the electronic spectrum of the 
carbon nanotubes, opening the possibility to tailor dif- 
ferent structures with the desired functionality. In the 
experiment reported in Ref. y, chains of Cqq molecules 
have been encapsulated in the interior of semiconducting 
nanotubes, showing that the measured spectra can be 
explained from the hybridization of the electronic states 
in the nanotubes with those in the fullerenes. In Ref. U, 
a different kind of peapod has been studied, formed by 
the encapsulation of sparse metallofullerenes Gd@Cg2 in 
semiconducting nanotubes. In that case, a striking obser- 
vation has been the spatial modulation of the nanotube 
bandgap, which is narrowed from 0.43 eV down to 0.17 
eV at the fullerene positions. 

Several electronic structure calculations have shown 
that some levels of the fullerene molecular orbitals may 
lie within the gap of the semiconducting nanotubes |g, Q- 
These results have been used to explain the narrowing of 
the bandgap observed experimentally by claiming that, 
at the fullerene sites, the measured gap should corre- 
spond to the distance from the top of the valence band 
to the lowest unoccupied molecular level within the nan- 
otube gap. This kind of effective reduction of the gap 
seems to be in agreement with the value measured in 
the experiments but, on the other hand, it leaves unac- 
counted a most significant property of its spatial modula- 
tion. Thus, while the effect of the unoccupied molecular 



orbital should decay over a distance of about 3A along 
the nanotube, the spatial modulation measured in Ref. 4 
is much less abrupt, with the narrowing of the gap taking 
place over an extension of about 2 or 3 ran. 

The distinctive feature of the peapods displaying the 
modulation of the bandgap lies in the charge of the en- 
capsulated fullerenes. For the nanotubes with arrays of 
Cqq molecules inside reported in Ref. 0, no sign of reduc- 
tion of the gap has been observed. On the other hand, it 
is natural to expect that the Coulomb interaction due to 
charged fullerenes may lead to significant effects in the 
electronic spectra. In that respect, it has been shown 
that, in the case of the peapods with the encapsulated 
metallofullerenes, there should be a charge transfer of 
about one electron from the latter to the surrounding 
nanotube. Thus, the metallofullerenes seem to have a 
net positive charge in the samples reported in Ref. 0, 
making pertinent the discussion of the Coulomb interac- 
tion. 

In the present paper, we show that the electrostatic in- 
teraction produced by charged fullerenes is able to mod- 
ify locally the electronic density of states of the peapods. 
In the cases where the interaction is felt as an attrac- 
tive potential by the electrons in the nanotube, localized 
electronic states are formed in the nanotubes around the 
position of the fullerenes. Their spatial extension as well 
as their number can be estimated by scaling arguments 
following from the particular form of the potential. For 
that purpose, we can rely on a one-dimensional (ID) low- 
energy model of the semiconducting nanotubes, which 
will allow also the determination of the bound states 
formed by the Coulomb potential within the gap. 

We illustrate our ID low-energy description in the case 
of the zig-zag nanotubes. These have a unit cell with 
length 3a, in terms of the C — C distance a . For each 
transverse array of atoms n = 1,2, ... N at a given nan- 
otube section x, it is convenient to introduce the Fourier 
transform of the electron operators c(x, n) 

c(x,n) ~ J2d p (x)e l2mlp / N (1) 
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where the index p labels the different subbands p — 
0, ±1, . . . ± N/2. Within the tight-binding approach, the 
hamiltonian for the zig-zag nanotube is 



H t b = — 1|| ^ dp (3la)dp(3la + a) 
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—t± z p dp(3la + a)d p (3la + 3a/2) 
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d p( 3la + 3a/2)d p (3/a + 5a/2) 
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— to. ^ z* p d+ (3/a + 5a/2)d p (3(/ + l)a) (2) 

= 1 + exp(i2np/ N) and 2 € Z. Henceforth 
we will consider that the transfer integrals for the two 
different bond orientations can be taken approximately 
equal, t± = t|| = t. Within this description, each band 
is obtained from the diagonalization of a one dimen- 
sional chain with two orbitals per unit cell, with gap 
2A p = 2t|l — 2 cos(2?rp/iV)|. This simple one-particle 
description accounts for the semiconducting character of 
the zig-zag nanotubes when N is not a multiple of 3. 

In the semiconducting nanotube, we can write N = 
3N' ± 1. The subbands closest to the Fermi level have 
p ~ N' , and a dispersion typical of a Dirac spinor 
with Fermi velocity v F = 3ta/2 and a mass parameter 
to = A p , e(k) w ±y / v F k 2 + to 2 Q. In the following we 
assume that the fullerene potential has radial symme- 
try, and the subband index p remains a good quantum 
number. Then, the calculation of the perturbed energy 
levels can be factorized into a set of equivalent problems, 
one for each subband of the nanotube. For the sake of 
describing the many-body effects in the interaction with 
the fullerene charge Q z , it is most convenient to write the 
model for the low-energy electronic states in the contin- 
uum, in terms of the Dirac spinor field with components 
^ l and ^>r. We will assume that the fullerene charge has 
a finite extension as seen along the longitudinal dimen- 
sion of the nanotube, over a distance R of the order of the 
fullerene radius. The full hamiltonian of the interacting 
theory becomes then 

H = J dx{-iv F ^+d x ^ L +iv F ^+d x ^ R 

+m*J*ij + 771^^ L ) 

e 2 f 1 
+— J dxdy(p L (x) + pr(x)) _ (p L (y) + Pr{v)) 



e 2 / dxdy 



1 



1 D -(x-zf/R 2 



x i \{Pl{v) +Pr{v)) 

\x-y\ 



(3) 



where pl(x), pr{x) are the respective electron density op- 
erators for the two spinor components 8J. 

The effect of the fullerene charge in the electronic spec- 
trum can be understood by applying the semiclassical 



approximation to the electron wavefunctions. In this ap- 
proach, we assume that the only significant effect of the 
electron-electron interaction within the nanotube is to 
induce some charge distribution screening the net charge 
of the fullerene. The screening is not complete, as it is 
produced by a system with a gap at the Fermi level, and 
it can be encoded at a later stage in the form of a suit- 
able dielectric function. At this point, we concentrate 
therefore on the resolution of the eigenvalue problem 

(-iv F d x +V(x- z))-® L + m$>R = £* L (4) 
(iv F d x + V(x - z))S> R + = ef R (5) 

The system (@J)-© can be reduced to an equation for 
any of the two components of the Dirac spinor, and the 
WKB solutions can be obtained in the usual fashion. 
In the semiclassical approximation, the \& r component 
reads for instance 



= ±i f x dxdxV(x-z)/^J(s-V(x-z)) 2 -m 2 
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±ij* dx(l/v F )y/(s-V(x-z)) 2 -m 2 
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The turning points are given by the roots of the equa- 



tion (e - V(r )) 2 = 



When such points exist, the 



semiclassical approximation gives rise to the quantiza- 
tion condition 



dx\/ (e — V(x)) 2 — Tn 2 = 
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(7) 



When the potential V(x) has the long-distance be- 
havior V(x) ~ — 1/| x|, it can be shown from Eq. (0 
that there is an infinite series of bound states with en- 
ergy below and approaching the value of to. If we define 
to — e = e, the quantization condition translates into 

r o / i 

dxy/(-V(x) - e)(2m - V(x) - e) = v f tt In + - 

r V 2 

(8) 

As the parameter e goes to zero, we observe that the 
integral at the left-hand-side of the equation scales as 
1 / y/e, reflecting the divergent range of integration. Thus, 
for large values of n, we find the quantization rule 
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The above arguments convey that the Coulomb po- 
tential of the fullerene charge has to induce localized 
electronic states within the nanotube gap. We show 
next that this conclusion holds after incorporating the 
screening of the potential by the electronic charge in the 
nanotube. Our continuum model allows actually an ex- 
act treatment of such effect by means of bosonization 
methods 9]. Following these, the electron density for a 
given subband near the Fermi level can be mapped into 
a boson field <j>(x) according to the relation -^d^(x) = 
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Pl{x) + Pr(x). In terms of the boson fields, the hamilto- 
nian @ can be expressed in the form 



H = 
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where II(x) is the momentum conjugate of the field $(cc) 
and a is a short-distance cutoff needed to regulate the 
product of fields at the same point. 

From the representation (|10|l , we observe that the pres- 
ence of the charge Q z induces a nonvanishing average 
value of the field d x $(x). As we are interested in low- 
energy screening properties, we can expand the cos(<i>(a;)) 
dependence to deal with a hamiltonian quadratic in the 
boson fields. The screening charge distribution can be 
obtained then by shifting the field $(x) by a suitable 
function, <§>{x) — f(x) = &(x), so that the redefined den- 
sity has a vanishing average, (d x $>(x)) = 0. The ap- 
proach allows an straightforward generalization to the 
case of Nf different electron flavors, accounting for spin 
and subband degeneracies. Thus, the expression for the 
screening electron density induced in the nanotube is in 
general 
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where the logarithmic dependences come from the rep- 
resentation of the bare Coulomb potential in momentum 
space [lfj. The shape of the screening charge distribution 
has been plotted in Fig. ^for m/t — 0.1 and R = 3a, 
taking Nf — 4 and e 2 /vf ~ 2.7 . We observe that in 
general there is only a partial screening of the fullerene 
charge Q z , due to the absence of gapless excitations in 
the semiconducting nanotube. The results are not qual- 
itatively changed if other subbands with higher gaps are 
also included. 

The dressed Coulomb potential Vd(x — z) felt in the 
nanotube is created by the sum of the background charge 
Q z e and the screening charge — ep s (x). Thus we have 



V d (x - z) 



dk 
2^ 



cos(k(x- z))e- R2k ^ 4 



(fc 2 + ^)21og(l + k /k) 



(12) 



We observe that the presence of the gap (accounted for 
by the m parameter) precludes an effective reduction of 
the Coulomb interaction at small momenta. The dressed 
potential Vd(x) (for the same choice of parameters made 
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FIG. 1: Plot of the model fullerene charge distribution for 
Q z — 1 (dashed line) and of the screening charge distribution 
induced in the carbon nanotube (full line). The variable x 
is measured in units of the C'-C distance in the nanotube. 
Inset: Log- log plot of the bare Coulomb potential (dashed 
line) and of the dressed potential obtained after considering 
the screening charge in the nanotube (full line). 



above) has been represented in the inset of Fig. ^ where 
it is also compared with the bare Coulomb potential. 

From the above analysis, it is clear that the dressed in- 
teraction keeps the long-range character of the Coulomb 
potential. Thus, in the more complete picture accounting 
for the screening effects within the nanotube, the result of 
having localized states due to the presence of the fullerene 
charge remains unmodified. The actual number of them 
depends on the value Q z of the net fullerene charge, as 
well as on the length of the carbon nanotube. The form of 
the spectrum within the gap can be obtained for a par- 
ticular system by diagonalizing the lattice hamiltonian 
made of the tight-binding term and the interaction 
term with the dressed potential 



Hi, 



V d {x l - z)d+(xi)d p (xi) 



(13) 



where Xi label the different sections of atoms in the nan- 
otube, in the same fashion as in Eq. 

The results of diagonalizing the hamiltonian H t b + Hi nt 
for a semiconducting (17,0) nanotube have been repre- 
sented in Fig. |2 The level of the deepest state within 
the gap depends on the net charge Q z , but there is always 
an accumulation of states near the conduction band edge. 
Using the continuum approximation, a state within the 
gap, at en ergy —A < e < A, is localized within a length 
I fa vp/\/ A 2 — e 2 . This estimate agrees well with the 
numerical calculations based on the tight-binding model. 
Hence, the localization length of the deepest levels within 
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FIG. 2: Diagram of the low-energy levels obtained from the 
diagonalization of the lattice hamiltonian for Q z — 0.5 (left) 
and Q z = 1.0 (right). The spectra are given in each case for 
a couple of nanotube segments with lengths L — 300a and 
1500a. The energies are in units of the transfer integral t. 
The dashed lines mark the border of the gap in the absence 
of the fullerene charge. 



the gap is comparable to the radius of the nanotube. 
Note that each subband will give rise to a set of localized 
states within its gap. 

The existence of localized states within the gap of the 
semiconducting nanotube has a significant impact on the 
shape of the local density of states of the system. This 
depends on the value of the net charge Q z and, for each 
particular experimental setup, on the amount of screen- 
ing due to external charges. For the sake of establishing 
a comparison with the spectra measured in Ref. |4], we 
have adopted a phenomenological approach by fixing the 
effective value of Q z so that the deepest localized state 
is placed at the center of the nanotube gap. The local 
density of states obtained from the diagonalization of the 
lattice hamiltonian H t t + Hint f° r a (17,0) nanotube is 
represented in Fig. [21 We observe the similarity with the 
local variation of the gap measured in Ref. around the 
position of each fullerene cluster. In particular, we notice 
that the narrowing of the gap takes place over a distance 
of a few nanometers, in agreement with the spatial mod- 
ulation of the gap observed in the experiment. 

Our results show that the Coulomb interaction be- 
tween the fullerene clusters and the nanotube has to be 
taken into account to achieve a consistent explanation of 



the spatial modulation of the gap in the peapods. In 
this respect, our approach is complementary to those 
studies that have focused on the hybridization of the 
molecular orbitals of the fullerenes with the nanotube 
states. The present analysis is consistent with the fact 
that the narrowing of the gap has been measured in the 
peapods with metallofullerenes, while no gap reduction 
has been observed in the samples with encapsulated Cqq 
molecules. Our investigation points at the possibility of 
having control of the gap reduction in the peapods, by 
varying the net fullerene charge Q z . This effect should 
be confronted in suitably prepared samples, for the sake 
of making progress in the local bandgap engineering of 
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FIG. 3: Intensity plot of the local density of states obtained 
from the diagonalization of the lattice hamiltonian for a (17,0) 
nanotube with Q z = 0.6. The variable x is measured in units 
of the C-C distance in the nanotube. 



nanotube devices. 
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